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ABSTRACT 


The  problem  of  optimal  frequency  separation 
involves  the  simultaneous  separation  of  several  vibra- 
tion modes  even  where  only  maximal  frequency  separation 
between  the  two  lowest  frequencies  is  desired.  Conven- 
tional optimization  procedures  which  maximize  a single 
function  cannot  deal  with  such  a problem.  Earlier  attempts 
to  treat  this  problem  have  therefore  encountered  diffi- 
culty in  solution.  New  optimization  methodology  is, 
therefore,  introduced  which  allows  the  simultaneous 
separation  of  several  modes.  This  new  procedure  is 
coupled  with  a powerful  new  general  optimization  algorithm. 
Test  results  indicate  that  this  new  combined  procedure 
is  capable  of  the  reliable  solution  of  the  optimal  frequency 
separation  problem. 


1. 


I.  Introduction 

Bronowickl  et  al.,  [1,2]  in  their  extension  of  earlier  studies  of 
optimal  synthesis  "T"  ring  stiffened  cylindrical  shells  under  hydrostatic 
pressure  [3,5]  introduce  the  problem  of  optimal  natural  frequency  separation. 
Their  pioneering  effort  uncovered  several  difficulties  associated  with  such 
problems . 

Ref.  [1]  in  its  discussion  of  their  type  III  problem  (separation  of 
the  two  lowest  frequencies  with  primarily  axial  motion)  notes  that  their 
optimal  search  terminates  at  a design  where  the  second  primarily  "axial" 
frequency  has  a radial  component  that  is  equal  to  the  axial  component.  In 
the  type  II  problem  (separation  of  the  two  lowest  frequencies)  their  search 
terminates  at  designs  where  the  second  and  third  frequencies  are  essentially 
equal.  The  authors  of  [1]  recognize  the  difficulty  in  obtaining  an  optimal 
solution  to  the  axial  frequency  separation  problem  which  they  attribute  to  the 
need  for  an  improved  definition  of  primarily  axial  frequencies. 

The  difficulty  with  these  problem  types  is  that  a conventional 
mathematical  programming  formulation  and  solution  procedure  is  not  effective 
at  the  points  described  above.  In  problem  type  II  at  points  where  the  second 
and  third  lowest  frequency  are  essentially  equal  it  is  necessary  to  separate 
the  first  and  third  as  well  as  the  first  and  second  frequencies  in  order  to 
continue  the  search.  In  problem  type  III  it  is  necessary  to  suppress  growth 
of  the  radial  component  associated  with  the  second  primariliy  axial  frequency 
in  addition  to  increasing  its  separation  from  the  lowest  axial  frequency  in 
order  to  continue  movement  toward  the  optimum. 

In  addition  to  the  above  difficulties  Ref.  [1]  utilizes  an  optimiza- 


tion procedure  which  is  apparently  unreliable  on  the  problems  formulated  therein. 


This  report  presents  new  methodology  for  treating  such  problems. 

It  applies  this  method  to  problem  type  II  of  Ref.  [1]  as  well  as  a similar 
problem  so  as  to  test  its  effectiveness  and  to  gain  additional  insight  into 
optimal  frequency  separation.  This  work  utilizes  a new  and  powerful  mathe- 
matical programming  procedure  [6],  which  appears  to  overcome  the  reliability 
problem  associated  with  that  used  in  Ref.  [1]  thereby  allowing  generation 
of  more  accurate  data  on  the  characteristics  of  designs  with  frequency 
separation . 

In  addition,  this  improved  optimization  procedure  is  adapted  to 
a minimum  weight  problem  formulation  similar  to  that  of  Ref.  [3]  and  prob- 
lem type  I of  Ref.  [1]  (minimum  weight  with  minimum  natural  frequency  con- 
straint) to  demonstrate  its  effectiveness  on  a six  variable  form  of  this 
problem.  Previous  attempts  to  solve  this  six  variable  form  had  proven  to 
be  insuccessful  [3,7]. 


3. 


II.  Problem  Formulation 

Mathematical  programming  procedures  treat  the  problem: 
Find  those  values  x^  of  the  variables  x^  that  result  in 


f(x^)  - min  f(x^)  i * 1,2, ...I 


subject  to  the  constraints 


gj (x±)  1 0 


j - 1,2. ..J 


£ u 

Xi  — Xi  — Xi 


These  procedures  start  from  some  initial  point  x°  and  by  some  strategy, 
usually  based  on  the  local  properties  the  functions  involved,  generate  a 


sequence  of  points  where 


f(x^+1)  < f(x*)  . 


In  problem  type  II  of  Ref.  [1] 

f (xf)  ■ Wj.  - u»2  (5) 

where  and  are  the  first  and  second  in  vacuo,  natural  frequencies  of  the 
shell  segment  shown  in  Figure  1.  The  variables  x^  are  the  skin,  stiffener 
web  and  stiffener  flange  thicknesses,  the  flange  width,  stiffener  spacing  and 
stiffener  web  height  respectively  as  given  in  Figure  1. 

Consider  the  problem  of  locating  an  improved  point  in  the  neighbor- 


hood of  a point  x'  where 


u2  ~ u3 


RING  AND  SKIN  DETAIL 


and  where  is  the  third  natural  frequency.  Each  frequency  where 
k - 1,2,3  - - « and 


^ u>2  “3  Z.  • • • • (6) 


is  a function  of  the  shell  design  parameters  such  as  its  radius,  length, 
material  properties,  applied  hydrostatic  pressure,  etc.,  the  problem  variables 

J r 

x^  (refer  to  Fig.  1)  and  the  number  of  circumferential  waves  n^  and  axial 

r r r 

half-waves  m^  associated  with  frequency  at  point  x^.  Thus,  for  a given 

th  t 

set  of  design  parameters  the  k natural  frequency  at  point  x^  is 


r _ , r r rx 

\ “(V  V V 


(7) 


a function  of  both  continuous  and  integer  variables.  Gradient  based  pro- 
cedures such  as  those  used  in  [1]  determine  a direction  for  function  im- 

r r it 

provement  from  the  derivates  of  and  with  respect  to  x^.  An  attempt 

to  move  in  a direction  increasing  the  separation  between  the  modes  associated 

with  these  frequencies  where  will  often  reduce  the  separation  be- 

r r 

tween  the  modes  associated  with  10^  and  These  modes  will  then  produce 

the  two  lowest  frequencies  after  a move  to  point  x^  where  now  the  separa- 
tion between  these  frequencies  will  be  lower  than  at  point  x£.  Conventional 
gradient  based  procedures  are  therefore  likely  to  fail  at  such  points. 

A procedure  is  needed  at  such  points  that  will  separate  from 
while  simultaneously  separating  and  Conventional  gradient  based 

procedures  cannot  provide  this  simultaneous  separation. 

Direct  search  optimization  procedures  avoid  this  difficutly  since 
they  do  not  involve  derivatives.  These  procedures  unfortunately  are  not 


highly  reliable  [8]  and  the  best  of  them  is  apparently  not  capable  of 
treating  the  difficult  six  variable  formulation  of  this  shell  design  problem 
[3,9].  On  the  other  hand  a combination  direct  search  and  gradient  based 
procedure,  the  Direct  Search-Feasible  Direction  (DSFD)  algorithm  [10] 
which  proved  reliable  on  a series  of  10  test  problems  of  Ref.  [8]  also 
proved  to  be  capable  of  treating  this  six  variable  shell  problem  [11],  The 
DSFD  thus  appears  to  be  adaptable  to  the  frequency  separation  problem. 

This  procedure  couples  an  efficient  direct  search  scheme  with  a gradient 
based  direction  finding  procedure  used  at  points  of  direct  search  failure. 

A refined  version  of  this  procedure  is  used  here  [6]. 

As  the  direct  search  avoids  the  use  of  derivatives  the  objective 
function  given  by  equation  (5)  can  be  used  without  difficulty  at  all  points 
where  direct  search  is  utilized.  At  points  of  direct  search  failure  however 
the  search  is  restarted  on  the  basis  of  local  derivative  information.  Fortun- 
ately the  direction  finding  algorithm  used  can  easily  be  modified  to  handle 
the  need  to  separate  and  as  well  as  and 

The  direction  finding  procedure  of  Zontendijk  [12]  is  utilized  to 
restart  the  basic  direct  search.  A direction  s^  in  which  an  improved  point 
may  be  found  is  normally  determined  from  the  solution  of  the  linear  pro- 
gramming problem: 

Given  the  set  x^,  find  the  set  s^  that  results  in  a 


for  which 


max  o 


o > 0 


(8) 

(9) 


(8l)T  V f(xt>  + a < 0 


(10) 


7. 


-(s±)  v gj.CXj)  + Wj  a 1 0 


-1  < 3^  - <_  1 


-1  < s,  < 0 i £ - I~ 
— i — a 


0 < s . < 1 i £ + I 
— i — a 


T 

Here  (s^)  indicates  the  transpose  of  vector  s^,  is  a weighting  parameter, 

the  set  J contains  the  active  constraints  where 
a 

gj  (x±)  < e1  (15) 

e^  is  small  arbitrary  positive  constraint  defining  "activity",  and  Ia>  I* 
constitute  the  active  upper  and  lower  regional  constraints,  respectively. 


u 

Xi  " Xi 


< e i £ I+ 


x . - x < e i £ I 
x i a 


Zoutendijk  [12]  shows  that  if  the  solution  s^  is  a null  vector, 
then  the  point  is  a local  optimum  and  if  not  then  the  direction  of  s^  is 
the  best  feasible  direction. 


Equations  (8-14)  constitute  a linq^k  programming  problem  with  the 
variables  s^  and  a.  Equation  (8)  is  the  objective  function  and  the  remaining 
equations  of  the  constraints.  The  solution  s^  can  be  obtained  reliably 
and  efficiently  using  any  suitable  linear  programming  method. 

At  points  where 

“3  “ “2  — e3  (17) 


8. 


and  e3  is  an  arbitrary  small  constant  indicating  frequency  similarity 
equation  (10)  is  replaced  by  the  set 


(s^1  VfJ(x1)+olO 


(3ir  Vf2(xi)  + a < 0 


(18) 


where  is  given  by  eq  (5)  and 


f2  " “l  “ “3 


(19) 


It  may  be  seen  that  the  solution  to  the  modified  direction  finding  problem 
will  yield  a direction  that  will  separate  cu^  and  u>3  as  well  as  separating 


and 


In  the  event  " a'3  * the  same  modification  procedure  may  be 


used.  Thus  if 


‘ “2  < e3 


(20) 


replace  equation  (10)  with  the  set 


where 


(si)  7 fl(xi}  + a < 0 


(s±)T  * f2(XjL)  + a < 0 ? 


(st)  V f3(xt)  + a < 0 J 


f3  " “l  ‘ “4 


(21) 


(22) 


Similar  reasoning  may  be  adapted  to  treat  problem  type  III  of 


Ref.  [1], 


9. 


III.  The  Shell  Design  Problem 

An  optimal  shell  design  capability  was  developed  along  the  lines  of 
earlier  procedures  [3,4]  to  allow  a study  of  the  aforementioned  optimal 
frequency  separation  method  and  to  provide  greater  design  flexibility  than 
available  in  earlier  programs. 

Four  problem  types  are  treated: 

Type  0 Problem 

The  objective  function  is 

f(x±)  = WD  (23) 

where  the  design  variables  are  as  given  in  Fig.  1 and  WQ  is  the  weight/dis- 
placement ratio  of  the  shell  segment  excluding  the  bulkheads  [3]. 

The  behavior  constraints  used  control; 

g^  = gross  buckling 

~ shell  (inter-ring)  buckling 
gg  = shell  yielding 

= stiffener  yielding 
gg  = stiffener  flange  buckling 
gg  = maximum  flange  thickness 
g -j  = minimum  flange  width 

gg  = minimum  internal  or  maximum  external  radius 

811  = We^  buckHn8 

constraints  gg  and  g^Q  are  not  used  with  the  type  0 problem. 

The  equations  used  for  gg,  g^,,  g&,  g7  and  gg  are  taken  from  [3]. 

Those  used  for  gg  and  g^g  are  essentially  eqns  (11)  of  Ref.  [2] . Reference 


A 


10. 


[2]  uses  the  same  equations  as  [3]  for  the  shell  and  stiffener  yield  con- 
straints and  except  that  [2]  ignores  the  effect  of  eccentricity  (out- 
of-roundness)  on  stiffener  stress. 

The  basic  behavior  prediction  equations  for  constraints  g^  and  g^ 
are  adapted  from  Ref.  [5]  which  uses  a procedure  described  in  [13].  Gross 
buckling  control  is  achieved  by  stating  that 


P - FS.P 
cr  1 

: k 

Per 


> 0 


(24) 


where  FS^  is  the  factor  of  safety  for  gross  buckling  p is  the  hydrostatic 


pressure  and 


p «*  min  p (n,m) 
rcr  cr  ' 


(25) 


where  n is  the  number  of  circumferential  waves  and  m the  number  of  axial 

half-waves.  The  equation  for  Pcr  and  the  method  for  finding  the  minimum 

buckling  load  p are  given  in  Ref.  [5].  This  equation  is  based  on  the 
cr 

Donnell  shell  theory  and  a smeared  stiffener  orthotropic  analysis. 

The  equation  for  the  shell  buckling  constraint  is  similar  to  that 
for  gross  buckling  except  that  x^,the  distance  between  stiffeners,  replaces 
L the  overall  shell  length  in  eqn  (2)  of  Ref.  [5]  and  the  stiffener  terms 

b 

are  omitted. 


Type  1 Problem 

This  problem  type  is  identical  to  type  0 except  a minimum  natural 
frequency  constraint  is  used  where 


w.  - w . 

1 min 


> 0 


(26) 


where  m ^ is  the  specified  minimum  frequency  and  the  lowest  natural 
frequency  of  the  structure.  The  quanity  is  found  from 

o)^  ■ min  w(n,m)  (27) 

where  m is  given  by  eqns  (31) —(33)  of  Ref.  [13]. 

The  n and  m producing  the  minimum  frequency  are  located  in  exactly 
the  same  fashion  as  for  the  minimum  buckling  load.  The  minimum  natural 
frequencies  of  both  the  entire  shell  and  the  segment  of  the  shell  between 
stiffeners  are  determined  and  the  lowest  value  is  used  in  eqn  (27). 

The  frequency  equations  are  derived  using  the  same  smeared  stiffen- 
er orthotropic  shell  approach  used  for  the  general  buckling  equation. 

The  effects  of  stiffener  placement  (interior  or  exterior)  and  torsional 
rigidity  are  considered  but  imperfection  sensitivity  is  ignored.  Experi- 
mental results  indicate  that  such  equations  produce  reasonably  accurate  re- 
sults for  hydrostatically  loaded  stiffened  shells  of  the  type  studied  here 
where  the  stiffener  spacing  is  small  compared  to  the  buckling  wave  length 
as  is  the  case  for  the  range  of  parameters  studied  in  this  report.  Imper- 
fections do  not  pl^v  a major  role  in  such  shells  [14,15]. 

It  should  be  noted  that  the  use  of  an  in  vacuo  model  for  the 
vibration  of  shells  submerged  in  water  may  produce  substantial  error  [7,16]. 

Its  use  here  however,  is  justified  in  light  of  the  objectives  of  this  work. 

These  objectives  are  to  present  and  evaluate  the  effectiveness  of  new 
methodology  for  the  design  of  shells  with  optimal  frequency  separation  and 
to  gain  insight  into  the  characteristics  of  such  shells. 

L * A 
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A comparison  of  the  values  of  the  inter-ring  buckling  pressures 
produced  by  equation  (8)  of  Ref.  [3]  and  the  procedure  used  here  shows  values 
within  2 Z for  the  designs  of  this  study.  Since  the  validity  of  the  equation 
used  in  [3]  has  been  demonstrated  by  experiment  [17]  the  model  used  here  for 
inter-ring  buckling  appears  to  be  suitable  even  though  it  assumes  simple 
support  of  the  ends  of  the  inter-ring  panel.  The  effect  of  this  assumption 
on  the  natural  frequency  values  for  such  short  shells  is  unknown.  The  use 
of  this  model  for  natural  frequency  prediction  in  this  preliminary  study, 
however,  seems  justified  in  light  of  the  rather  substantial  error  produced 
by  the  in  vacuo  assumption  and  the  objectives  of  the  study.  A consideration 
of  inter-ring  vibration  allows  an  opportunity  to  determine  if  this  mode 
appears  to  be  significant  for  the  types  of  problems  studied  here.  Further- 
more inclusion  of  this  vibration  mode  avoids  the  inconsistancy  in  Ref.  [1] 
wherein  both  an  inter-ring  buckling  constraint  (natural  frequency  is  zero 
where  such  a constraint  is  active)  and  a minimum  natural  frequency  constraint 
are  specified  [7]. 


Type  2 Problem 

Here  the  objective  function  is  given  by  equation  (5)  and  a maximum 
weight  constraint 


s10 


W - 
max  D 


> 0 


(28) 


is  used  in  addition  to  gQ,  where  W is  the  specified  maximum  allowable 

7 QlEX 

weight/displacement  ratio. 


This  is  the  optimal  frequency  separation  problem  (type  II)  of 


Type  3 Problem 


The  objective  in  this  problem  type  is  to  maximize  the  lowest 
natural  frequency.  Thus  the  objective  function  here  is 

f(x^  — — ou1  (29) 

This  problem  type  uses  the  same  constraints  as  type  2. 


A difficulty  similar  to  that  arising  in  problem  type  2 is  encount- 
ered here.  The  optimal  search  will  encounter  points  where  3 idj  3 

This  situation  is  treated  in  the  same  fashion  as  for  the  type  2 problem. 
Thus  where 

let  w-  - < e^  (30) 


(31) 


and  equations  (18)  are  used  in  place  of  (10).  Where 


let 

f3  * ' “3 

and  replace  eqn  (10)  by  eqns  (21) . 


(32) 


(33) 


Two  additional  options  are  available  to  the  user  of  this  capability. 
One  option  allows  the  designer  to  use  a five  variable  formulation.  Here 
the  value  of  the  web  thickness  is  calculated  so  as  to  just  satisfy  the  web 
buckling  constraint.  Thus  web  thickness  is  not  treated  as  an  independent 
variable  when  this  option  is  selected.  Since  the  web  buckling  constraint 
was  found  to  be  active  in  all  cases  studied  this  option  can  apparently  be 
used  to  reduce  computational  effort  while  still  producing  optimal  solutions. 


14  . 


Another  option  specifies  a four  variable  formulation.  Here  the 
number  of  stiffeners  to  be  used  is  specified  rather  than  treated  as  a 
variable  and  the  web  thickness  is  calculated  as  above . Since  it  was  found 
in  earlier  studies  [3,4]  that  there  is  usually  little  increase  in  weight 
associated  with  specifying  an  arbitrary  stiffener  spacing,  the  designer 
is  thus  for  practical  purposes  free  to  choose  a stiffener  spacing,  within 
relatively  broad  limits,  with  little  performance  penalty  for  problem  types 
0 and  1. 

The  four  variable  formulation  may  be  used  with  the  type  2 and  3 
problems  but  the  impact  on  performance  resulting  from  arbitrarily  selecting 
stiffener  spacing  has  not  been  studied. 


IV.  Results 

A computer  program  called  SBSHL7  coupling  the  DSFD  optimization 

procedure  and  problems  described  above  was  used  to  repeat  some  of  the  shell 

design  studies  of  Refs.  [1,3]  in  order  to  evaluate  the  performance  of 

SBSHL7  on  this  problem,  to  investigate  the  problem  of  optimal  frequency 

separation,  and  to  allow  a comparison  of  the  results  generated  with  those 

of  Refs.  [1,3].  All  studies  use  the  following  design  parameters  unless 

otherwise  indicated;  shell  midplane  radius  R * 5.029m  (198  in.);  shell 

length  Lg  - 15.09m  (594  in.)  shell  eccentricity  of  zero;  immersion  depth 

304.8m  (1000  ft);  specific  weight  of  immersion  fluid  (sea  water)  y =*  1.0256 

w 

3 3 3 

g/cm  (0.03741b/in  ),  specific  weight  of  shell  material  yg  ■ 7.733  g/cm 

(0.282  lb/in. ^)  Young's  modulus  20.68  X 10^N/cm^  (30  X lO^psi) , Poisson's 

2 

ratio  y ■ 0.30  and  allowable  yield  stress  of  41,360  N/cm  (60,000  psi)  for 
the  shell  and  stiffener  material.  Factors  of  safety  of  2 are  used  for  all 
buckling  constraints  and  a minimum  natural  frequency  of  zero  is  specified. 

All  shells  use  interior  stiffeners. 

Table  1 demonstrates  the  reliability  of  the  new  optimization  procedure 
on  the  type  0 problem.  The  starting  points,  in  mm,  used  for  this  study  were: 

(xj)x-  (0,0, 0,0, 0,0) 

(x°)2-  (12.7,12.7,12.7,127,254,127) 

(x°)3«  (25.4,25.4,25.4,254,508,254) 

(x°)4«  (38.1,38.1,38.1,381,762,381) . 


r 
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Table'  1.  Convergence  by  SBSHL7  to  an  Optimal 
Widely  Separated  Starting  Points  for  the  Type 


Configuration 
0 Problem. 


Ref.  [8] 

Starting 
Point  1 

Starting 
Point  2 

Starting 
Point  3 

Starting 
Point  4 

w 

V 

weight/ 

displacement  ratio 

0.1357 

0.10306 

0.10303 

0.10313 

0.10312 

V 

skin  thickness , mm  (in) 

26.901 

(1.0591) 

28.092 

(1.1060) 

28.128 

(1.074) 

28.082 

(1.1056) 

28.067 

(1.1050) 

V 

web  thickness,  mm  (in) 

24.226 

(0.9538 

6.7158 

(0.2644) 

6.6751 

(0.2628) 

6.7539 

(0.2659) 

6.7005 

(0.2638) 

x3. 

flange  thickness,  mm  (in) 

9.9772 

(0.3928) 

33.609 

(1.3232) 

27.277 

(1.0739) 

12.893 

(0.5076) 

15.133 

(0.5958) 

X4* 

flange  width,  mm  (in) 

125.308 

(4.933) 

25.293 

(0.9958) 

30.322 

(1.1938) 

67.774 

(2.6683) 

55.535 

(2.1864) 

V 

stiffener  spacing,  mm  (in) 

580.29 

(22.846) 

405.79 

(15.9763) 

407.67 

(16.050) 

412.47 

(16.239) 

406.30 

(15.996) 

V 

web  height,  mm  (in) 

463.07 

(17.168) 

292.68 

(11.5238) 

296.93 

(11.690) 

299.14 

(11.777) 

298.37 

(11.747) 

81> 

gross  buckling 

0.622 

0.000 

0.000 

0.000 

0.000 

82  * 

panel  buckling  ^ 

0.003 

0.440 

0.438 

0.423 

0.437 

8-j’ 

skin  yield 

0.004 

0.000 

0.000 

0.000 

0.000 

®4' 

stiffener  yield 

0.339 

0.088 

0.088 

0.089 

0.089 

85» 

flange  buckling 

0.003 

1.000 

0.993 

1.000 

0.989 

g6* 

web  buckling 

* 

0.041 

0.001 

0.013 

0.000 

*Web  buckling  controlled  by  setting  x^  ■ 18  x^. 


It  may  be  seen  that  a series  of  runs  from  widely  separated  starting 
points  converged  to  similar  designs  with  weights  that  are  within  0.15%  of  the 
lowest.  This  behavior  is  indicative  of  reliable  optimization  algorithm  per- 
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formance  and  the  absence  of  local  optima. 

The  designs  generated  by  SBSHL7  are  substantially  lighter  than  those 
given  in  Ref.  [3]  due  primarily  to  the  elimination  of  an  unrealistic  stiff- 
ener buckling  constraint  used  in  [3]  but  not  used  here  or  in  [1,5]. 

The  shell  buckling  equations  used  here  are  considerably  more  com- 
plex than  those  of  [3].  This  added  complexity  in  addition  to  the  increase 
in  problem  dimensionality  and  difficulty  resulting  from  the  uncoupling  of 
the  stiffener  variables  greatly  increases  the  computational  effort  required 
for  solution.  The  program  SBSHL7  typically  requires  about  100  times  more 
CPU  time  than  that  of  [3].  Thus,  the  added  sophistication  is  obtained  at 
substantially  increased  cost.  Still,  if  a reasonable  starting  point  is 
specified,  the  cost  of  solution  is  not  excessive  (normally  less  than  2 min 
CPU  time  on  an  IBM  370-168).  The  buckling  equations  used  in  [3]  can  un- 
fortunately lead  to  invalid  designs  for  certain  ranges  of  parameter  values 
[5]  since  the  interior  buckling  minimum  (m>l)  is  not  considered.  This  is 
also  true  of  the  procedure  of  Ref.  [1]  where  m>6  is  Ignored.  Thus,  this 
added  complexity  is  justified.  Fortunately  the  computational  cost  can 
apparently  be  drastically  reduced  by  use  of  a five  variable  formulation 
as  will  be  shown  below. 

Table  2 demonstrates  the  superiority  of  DSFD  over  the  SUM!  procedure 
used  in  [1]  which  failed  to  locate  the  minimum,  producing  designs  substantially 
heavier  than  those  presented  here.  The  shell  parameters  utilized  in  [1] 


r 
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Table  2.  Comparison  of  Designs  Developed  by  SBSHL7 
With  Those  of  Refs.  [1,2]  for  the  Type  1 Problem. 


Ref.  [1] 

Ref.  [2] 

SBSHL7 

Wp,  weight/ 

displacement  ratio 

0.13317 

0.11274 

0.10746 

x^,  skin  thickness,  mm  (in) 

30.754 

(1.2108) 

30.622 

(1.2056) 

29.522 

(1.1623) 

web  thickness,  mm  (in) 

9.563 

(0.3765) 

6.0274 

(0.2373) 

5.1054 

(0.2010) 

x^,  flange  thickness,  mm  (in) 

11.951 

(0.4705) 

7.8003 

(0.3071) 

4.6533 

(0.1832) 

x^,  flange  width,  mm  (in) 

448.66 

(17.664) 

263.22 

(10.363) 

150.80 

(5.937) 

x^,  stiffener  spacing,  mm  (in) 

853.49 

(33.602) 

766.32 

(30.170) 

424.51 

(16.713) 

Xg,  web  height,  mm  (in) 

497.56 

(19.589) 

279.91 

(11.020) 

239.04 

(9.411) 

g^,  gross  buckling 

Not  given 

Not  given 

0.234 

g^,  panel  buckling 

Not  given 

Not  given 

0.742 

g^,  skin  yield 

Not  given 

Not  given 

0.000 

g^,  stiffener  yield 

Not  given 

Not  given 

0.081 

g^ , flange  buckling 

Not  given 

Not  given 

0.039 

gg,  web  buckling 

Not  given 

Not  given 

0.000 

g_,  min.  nat.  freq.  H 
/ z 

Not  given 

Not  given 

0.000 

“1*  HZ 

28.12 

12.03 

12.00 

UI-,  H 

49.39 

22.30 

20.50 

► 


The  design  of  Ref.  [1]  is  obviously  not  optimal  since  an  optimal 
design  should  converge  to  the  minimum  frequency  constraint.  Neither  is 
the  design  a local  optimum.  This  design  was  used  as  a starting  point  for  a 
synthesis  run  using  SBSHL7.  The  search  immediately  located  a better  near- 
by design  and  moved  to  a design  essentially  identical  to^that  presented 
in  Table  2 for  SBSHL7.  Since  the  general  and  shell  buckling  constraints 
are  not  active  at  this  point,  the  constraint  equations  used  here  are  identical 
to  Ref.  [1].  Thus,  the  design  of  Ref.  [1]  does  not  appear  to  be  a local 
optimum. 


The  design  presented  in  Table  2 generated  by  SBSHL7  used  the  start- 
ing point  (x®)^ 


The  design  generated  in  Ref.  [2]  cannot  be  considered  as  representa- 
tive of  the  performance  of  the  SUMT  procedure  of  [1]  since  it  used  a near 
optimal  starting  point  supplied  by  Pappas  [7].  The  performance  of  SUMT 
in  [1,2]  coupled  with  the  evidence  developed  in  an  earlier  comparison  study 
[8],  strongly  suggests  that  SUMT  simply  cannot  cope  with  this  six  variable 
shell  design  problem. 


Problem  type  2 and  3 where  run  using  both  five  and  six  variable 
formulations.  Table  3 give  the  results  of  typical  problem  type  2 and  3 
runs.  A maximum  weight  displacement  ratio,  (WQ)  of  0.150  was  used  for  these 
runs.  The  quantities  n^  and  ng  are  the  number  of  objective  and  constraint 
function  evaluations  required  for  convergence.  Several  observations  can  be 


made  from  the  data  contained  in  this  table. 


Table  3. 

Optimal  Designs 

for  Problem  Type 

2 and  3 

Initial 

Point 

Type  2,  5 Var. 
Optimum 

Type  2,  6 Var. 
Optimum 

Type  3,  6 Var. 
Optimum 

Type  II 
Ref.  [1] 

WD 

0.155 

0.150 

0.150 

0.150 

0.135 

25.40 

38.46 

38.56 

35.00 

31.03 

X1 

(1.000) 

(1.514) 

(1.518) 

(1.378) 

(1.222) 

25.40 

13.97 

13.69 

11.81 

10.03 

X2 

(1.000) 

(0.550) 

(0.539) 

(0.465) 

(0.395) 

25.40 

15.29 

33.81 

23.04 

11.82 

x3 

(1.000) 

(0.602) 

(1.331) 

(0.907) 

(0.465) 

254.0 

246.4 

122.6 

541.0 

445.8 

X4 

(10.000) 

(9.702) 

(4.826) 

(21.30) 

(17.55) 

508.0 

1396 

1398 

1372 

859.9 

X5 

(20.000) 

(54.96) 

(55.04) 

(54.00) 

(33.85) 

254.0 

820.9 

804.9 

731.0 

526.3 

X6 

(10.000) 

(32.32) 

(31.69) 

(28.78) 

(20.72) 

81 

0.539 

0.789 

0.790 

0.749 

0.677 

g2 

0.472 

0.304 

0.307 

0.126 

0.342 

g3 

0.001 

0.114 

0.116 

0.001 

-0.002 

g4 

0.389 

0.467 

0.466 

0.521 

0.325 

g5 

0.947 

0.868 

0.994 

0.714 

0.619 

g6 

0.000 

0.602 

0.113 

0.341 

0.018 

g7 

10.000 

9.702 

4.826 

21.998 

17.551 

g8 

0.598 

0.544 

0.543 

0.553 

0.574 

g9 

0.534 

0.717 

0.717 

0.743 

0.673 

OQ 

M 

O 

-0.035 

0.001 

0.001 

0.000 

0.120 

gll 

0.967 

0.000 

0.000 

0.005 

-0.006 

0)1  (n1m1)Hz 

25.76(2,1) 

42.40(2,1) 

42.41(2,1) 

46.71(2,1) 

28.37(36.65)* 

(n2m2)Hz 

33.04(3,1) 

72.84(1,1) 

72.90(3,1) 

46.87(14,1)* 

51.96(58.52)** 

u>^(n^m^)Hz 

44.29(3,2) 

72.84(3,1) 

72.90(1,1) 

47.50(15,1)* 

51.96 

<Vn4m4)Hz 

47.82(1,1) 

72.88(14,1)* 

73.19(14,1)* 

49.83(13,1)* 

- 

(id2-u»1)Hz 

7.28 

30.44 

30.49 

0.16 

23.59(21.86)** 

1 

2963 

16.779 

22,572 

- 

n 

g 

11 

2118 

34,371 

50,193 

- 

*Mode  associated  with  shell  (inter-ring)  vibration. 

**Frequency  using  Donnell  shell  theory  employed  in  this  report . 
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The  computational  effort  required  to  achieve  convergence  using  the 
five  variable  formulation  is  substantially  less  than  that  using  the  six 
variable  form.  Since  all  studies  of  all  problem  types  indicates  that  the 
web  buckling  constraint  is  always  active  at  the  optimum  there  seems  to  be 
little  reason  to  use  the  six  variable  formulation  for  design  purposes.  The 
principal  advantage  of  this  form  is  that  it  provides  a costly  check  of  the 
assumption,  used  in  the  five  variable  form,  that  web  buckling  is  active  at 
the  optimum. 

An  examination  of  the  optimal  design  for  the  type  2 problem  reveals 
several  interesting  characteristics.  Not  only  is  essentially  equal  to 
1^2  as  found  in  Refs.  [1,2]  but  is  also  essentially  equal  to  Further- 

more, is  associated  with  shell  panel  (inter-ring)  vibration.  The  design 
is  characterized  by  relatively  large  frame  spacing  and  large,  deep  framing 
members.  Thus,  the  optimal  design  uses  the  largest  frame  members  and 
spacing  possible  without  inducing  a shell  panel  vibration  mode  lower  than 
the  second  frequency  associated  with  the  gross  vibration,  and  without  vio- 
lating the  maximum  weight  constraint.  It  is  apparent  therefore  that  the 
shell  panel  vibration  should  be  considered  in  optimal  frequency  separation 
problems . 

Since  the  designs  are  controlled  by  the  maximum  weight  constraint 
it  appears  that  optimal  separation  is  achieved  by  paying  a penalty  in  weight. 

The  characteristics  of  the  design  with  the  maximum  lowest  frequency 
is  rather  similar  to  that  with  the  largest  frequency  separation.  The  former 
has  a somwhat  higher  lowest  frequency  but  no  significant  separation  of  this 
frequency  from  the  second  or  third  lowest  frequencies.  Here  again,  shell 
panel  vibration  controls  the  design. 


Comparing  Che  designs  in  table  3 it  appears  that  rather  similar 
designs  may  behave  quite  differently  with  respect  to  frequency  separation. 

This  raises  a serious  question  with  regard  to  the  validity  of  frequency 
separation  data.  Orthotropic  shell  theory  is  only  an  approximation  to  actual 
behavior.  This  and  the  earlier  study  of  Refs.  [1,2]  furthermore  use  in 
vacuo  frequencies  to  study  the  characteristics  of  submerged  shells.  Con- 
sidering the  inaccuracy  of  these  assumptions  and  the  sensitivity  of  the 
frequency  separation  results  to  relatively  small  design  changes  the  data 
for  the  type  2 problem  appear  to  have  little  meaning  except  to  indicate 
possible  general  characteristics  of  such  designs. 

In  addition  to  the  search  starting  point  cited  in  Table  3 were 
made  using  the  six  variable  formulation  for  problem  types  2 and  3 from  the 
three  additional  starting  points  used  in  the  type  0 study.  These  also 
converged  to  designs  similar  to  those  of  table  3 thus  the  starting  point 
sensitivity  noted  in  Ref.  [1]  on  their  problem  type  2 was  not  apparent  here. 
These  results  indicate  that  the  optimization  procedure  used  here  is  apparently 
capable  of  locating  an  optimum  design  with  reasonable  reliability  for  all 
problem  types  studies . 

It  should  be  noted  that  the  frequency  separation  obtained  for  the 
type  2 problem  (u^  - - 30.48  Hz)  is  substantially  greater  than  obtained 

in  Ref.  [1]  (u^  - - 23.59  Hz)  and  that  Ref.  [1]  failed  to  converge  to 

the  maximum  weight  constraint.  It  appears  therefore  that  the  SUMT  procedure 
used  in  [1]  also  failed  to  locate  the  optimum  design  in  the  type  2 problem 
of  Ref.  [1].  The  difference  in  frequency  separation  results  do  not  appear 
to  be  due  to  the  difference  in  shell  theory  used  in  [1]  (Flugge-Lure-Byme) 
with  that  used  here  (Donnel  type).  The  point  cited  as  optimal  in  [1]  produced 
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. 

a separation  of  21.86  Hz  using  the  Donnel  theory  utilized  here,  a value 
somewhat  lower  than  that  of  [1],  but  still  in  good  agreement.  This  lower 
value  further  suggests  that  the  frequency  separation  at  the  optimal  point 
claimed  here,  would  in  fact  be  somewhat  larger  using  the  frequency  pre- 
diction procedure  of  [1], 
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IV.  Conclusion 

Design  for  optimal  frequency  separation  may  require  treatment  of 
the  case  where  several  frequencies  must  be  separated  simultaneously.  The 
methodology  presented  here  appears  effective  in  treating  such  problems. 

These  studies  indicate  that  the  simultaneous  separation  procedure  and  DSFD 
optimization  algorithm  utilized  here  are  capable  of  reliably  treating  the 
optimal  frequency  separation  problem. 

The  design  studies  show  that  shell  panel  vibration  modes  may  be 
active  at  the  optimum.  Thus,  a consideration  these  modes  should  be  included 
in  the  optimal  frequency  separation  problem.  The  similarity  of  designs  with 
optimal  and  negligable  separation  of  the  two  lowest  natural  frequencies  as 
predicted  by  the  orthotropic  in  vacuo  model  used  here  ana  elsewhere  indicates 
that  more  accurate  behavior  prediction  models  are  needed  to  produce  meaning- 
ful design  results. 
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